Methods for Selecting Medications

ABSTRACT

The invention provides computer assisted methods for selecting an antidepressant medication for a patient. The methods utilize a combination genotype comprising at least three genes to guide medication selection.

TECHNICAL FIELD

This invention relates to methods for selecting a medication (e.g., psychotropic medication) for a patient, and more particularly to selecting a patient's medication based on the genotype of genes encoding drug-metabolizing enzymes and genes encoding products involved in, for example, neurotransmission.

SUMMARY

The invention is based on the identification of a set of genes with polymorphisms that are associated with psychiatric diagnosis and pharmacological response to a medication. As a result, methods of the invention allow the genotype of a patient to be determined and, based on the genotype, a suitable medication to be selected for the patient. Methods of the invention allow the output of multiple genotypic assessments to be integrated, providing important and improved clinical information on which to select and dose medications. Thus, the methods of the invention provide a rational method for the identification of a medication that will result in an optimal response in the patient.

In one aspect, the invention features a method for selecting a psychotropic medication (e.g., an antidepressant, an antipsychotic agent, an antianxiety-sedative agent, or a mood-stabilizing agent) for a patient. The method includes providing the patient's genotype for a panel of genes, wherein the panel includes at least three cytochrome P450 genes and a serotonin transporter gene and selecting the psychotropic medication based on the genotype. The at least three cytochrome P450 genes can encode CYP2D6, 1A2, and 2C19 and providing the patient's genotype can include determining if the patient has the CYP1A2*1A or 1A2*3 allele, the CYP2C19*1A, 2C19*1B, or 2C19*2A allele, and the CYP2D6*1A, 2D6*2, 2D6*2N, 2D6*3, 2D6*4, 2D6*5, 2D6*6, 2D6*7, 2D6*8, 2D6*10, 2D6*12, or 2D6*17 allele.

The panel can include 4 cytochrome P450 genes (e.g., genes encoding CYP2D6, 2C19, 3A4, and 1A2). Providing the patient's genotype can include determining if said patient contains the CYP2D6*2, 2D6*3, 2D6*4, 2D6*10, or 2D6*17 allele, the CYP2C19*2A, 2C19*2B, 2C19*3, 2C19*4, 2C19*5A, 2C19*5B, 2C19*6, 2C19*7, or 2C19*8 allele, the CYP3A4*1B, 3A4*2, 3A4*5, 3A4*6, 3A4*12, 3A4*13, 3A4*15A, 3A4*17, 3A4*18A allele, and the CYP1A2*1F allele.

The panel can include at least 5 cytochrome P450 genes (e.g., CYP1A1, 1A2, 2D6, 2C19, and 3A4). Providing the patient's genotype can include determining if the patient has the CYP1A1*1A, 1A1*2, 1A1*3, or 1A1*4 allele, the CYP1A2*1A or 1A2*3 allele, the CYP2C19*1A, 2C19*1B, or 2C19*2A allele, the CYP2D6*1A, 2D6*2, 2D6*2N, 2D6*3, 2D6*4, 2D6*5, 2D6*6, 2D6*7, 2D6*8, 2D6*10, 2D6*12, 2D6*17, or 2D6*35 allele, and the CYP3A4*1A or 3A4*1B allele.

The panel further can include a plurality of dopamine receptor genes (e.g., dopamine receptor genes encoding dopamine receptors D1, D2, D3, D4, D5, and D6), a plurality of serotonin receptor genes (e.g., serotonin receptor genes encoding serotonin receptors 1A, 1B, 1D, 2A, or 2C), a tryptophan hydroxylase gene, and/or a catechol-O-methyl transferase gene. Serotonin receptor gene 2A is particularly useful.

In some embodiments, the panel includes at least 10 cytochrome P450 genes (e.g., genes encoding CYP1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, and 3A5). Determining the patient's genotype can include determining if the patient has the CYP1A1*1A, 1A1*2, 1A1*3, or 1A1*4 allele; the 1A2*1A or 1A2*3 allele; the CYP1B1*1, 1B1*2, 1B1*3, 1B1*4, 1B1*11, 1B1*14, 1B1*18, 1B1*19, 1B1*20, or 1B1*25 allele; the CYP2A6*1A, 2A6*1B, 2A6*2, or 2A6*5 allele; the CYP2B6*1, 2B6*2, 2B6*3, 2B6*4, 2B6*5, 2B6*6, or 286*7 allele; the CYP2C8*1A, 2C8*1B, 2C8*1C, 2C8*2, 2C8*3, or 2C8*4 allele; the CYP2C9*1, 2C9*2, 2C9*3, or 2C9*5 allele; the CYP2C18*m1 or 2C18*m2 allele; the CYP2C19*1A, 2C19*1B, or 2C19*2A allele; the CYP2D6*1A, 2D6*2, 2D6*2N, 2D6*3, 2D6*4, 2D6*5, 2D6*6, 2D6*7, 2D6*8, 2D6*10, 2D6*12, 2D6*17, or 2D6*35 allele; the CYP2E1*1A, 2E1*1C, 2E1*1D, 2E1*2, 2E1*4, 2E1*5, or 2E1*7 allele; the CYP3A4*1A or 3A4*1B allele; and the CYP3A5*1A, 3A5*3, 3A5*5, or 3A5*6 allele.

In some embodiments, the panel of genes contains the CYP2D6, 2C19, 3A4, and 1A2 genes, the serotonin receptor gene, and the serotonin receptor 2A gene. Selecting the psychotropic medication can include correlating genotype of the cytochrome P450 genes with capacity of each cytochrome P450 enzyme encoded by each cytochrome P450 gene to metabolize the psychotropic medication and correlating genotype of the serotonin transporter gene and the serotonin receptor gene with ability of the patient to respond to the psychotropic medication.

In another aspect, the invention features an article of manufacture that includes a substrate. The substrate includes a plurality of discrete regions, wherein each region includes a different population of nucleic acid molecules, wherein the different populations of nucleic acid molecules independently include nucleic acid molecules for detecting the CYP1A1*1A, 1A1*2, 1A1*3, and 1A1*4 alleles; the CYP1A2*1A and 1A2*3 alleles; the CYP2C19*1A, 2C19*1B, and 2C19*2A alleles; the CYP2D6*1A, 2D6*2, 2D6*2N, 2D6*3, 2D6*4, 2D6*5, 2D6*6, 2D6*7, 2D6*8, 2D6*10, 2D6*12, 2D6*17, and 2D6*35 alleles; the CYP3A4*1A and 3A4*1B alleles; the DAT1 40 bp VNTR and 10 repeat alleles; and the 5HTT promoter repeat and exon 2 variable repeat alleles. The population of nucleic acid molecules further can include nucleic acid molecules for detecting the CYP1B1*1, 1B1*2, 1B1*3, 1B1*4, 1B1*11, 1B1*14, 1B1*18, 1B1*19, 1B1*20, and 1B1*25 alleles, the CYP2A6*1A, 2A6*1B, 2A6*2, and 2A6*5 alleles, the CYP2B6*1, 2B6*2, 2B6*3, 2B6*4, 2B6*5, 2B6*6, and 2B6*7 alleles, the CYP2C8*1A, 2C8*1B, 2C8*1C, 208*2, 2C8*3, and 2C8*4 alleles, the CYP2C9*1, 2C9*2, 2C9*3, and 2C9*5 alleles, the CYP2C18*m1 and 2C18*m2 alleles, the CYP2C19*1A, 2C19*1B, and 2C19*2A alleles, the CYP2E1*1A, 2E1*1C, 2E1*1D, 2E1*2, 2E1*4, 2E1*5, and 2E1*7 alleles, and the CYP3A5*1A, 3A5*3, 3A5*5, and 3A5*6 alleles. The population further can include nucleic acid molecules for detecting the TPH A218C, A779C, G5806T, A6526G, and (CT)m(CA)n(CT)p alleles.

The invention also features a method of building a database for use in selecting a medication for a patient. The method includes receiving, in a computer system, a plurality of genotypes for a panel of genes; receiving a plurality of medication profiles specified based on the genotypes; and storing the plurality of genotypes and the medication profiles such that each medication profile is associated with one of the genotypes. The at least one medication profile can identify a medication and the medication can be placed in one of multiple categories included in the medication profile. Such categories can be selected from the group consisting of: medications that are safe to use, medications that should be used with caution, medications that should be closely monitored when used, medications that should be avoided, and combinations thereof. The medication profile can identify a universe of possible medications for the patient's genotype.

In another aspect, the invention features a computer program product containing executable instructions that when executed cause a processor to perform operations. The operations can include: receive a plurality of genotypes for a panel of genes; receive a plurality of medication profiles specified based on the genotypes; and store the genotypes and the medication profiles such that each medication profile is associated with one of the genotypes.

The invention also features a method of selecting a medication for a patient. The method includes receiving, in a computer system, a patient's genotype for a panel of genes; identifying, in a database comprising a plurality of medication profiles associated with genotypes, a medication profile that is associated with the patient's genotype; and outputting the identified medication profile in response to receiving the patient's genotype. A user can enter the patient's genotype in the computer system or the patient's genotype can be received directly from equipment used in determining the patient's genotype.

The medication profile can include a ranking of several medications, e.g., based on specific co-factors. The method can include adjusting the ranking before outputting the identified medication profile (e.g., based on receiving a genotypic polymorphism carried by the patient or based on receiving a clinical response relating to the patient). The clinical response can be by a family member of the patient.

In yet another aspect, the invention features a computer program product containing executable instructions that when executed cause a processor to perform operations that include receive a patient's genotype for a panel of genes; identify, in a database including a plurality of medication profiles associated with genotypes, a medication profile that is associated with the patient's genotype; and output the identified medication profile in response to receiving the patient's genotype.

Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although methods and materials similar or equivalent to those described herein can be used to practice the invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.

Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.

DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of a computer system 100, according to one embodiment.

FIG. 2 is a flow chart of a method 200 for building a database for use in selecting a medication for a patient.

FIG. 3 is a flow chart of a method 300 for selecting a medication for a patient.

DETAILED DESCRIPTION

In general, the invention features a method for selecting a medication (e.g., psychotropic medication) for a patient based on the genotype of genes that are useful for medication selection. Genes to be genotyped typically encode products that influence the metabolism of a medication or that are associated with a better treatment response. An algorithm can be used that initially excludes medications based on problematic genotypes of drug metabolizing enzyme genes (e.g., cytochrome P450 genes). A second stage of such an algorithm can begin with a list of potentially appropriate medications (e.g., antidepressants) and result in further classification of optimal treatments based on the genotype of target genes. For example, for antidepressants, further classification of the potential antidepressants can be based on the genotype of the serotonin transporter (5HTTR) and serotonin receptor 2A (HTR2A).

Psychotropic medications include antipsychotics or neuroleptics, antianxiety-sedative agents, antidepressants or mood elevating agents, and mood-stabilizing drugs such as lithium salts or valproic acid. Non-limiting examples of antipsychotics include the tricyclic phenothiazines (e.g., chlorpromazine, triflupromazine, thioridazine, and mesoridazine, fluphenazine, and trifluoperazine), thioxanthenes (e.g., chlorprothixene, clopenthixol, flupenthixol, piflutixol, and thiothixeue), and dibenzepines (e.g., loxapine, clozapine, clothiapine, metiapine, zotapine, ICI-204,636, fluperlapine, and olanzapine); butyrophenones such as haloperidol, diphenylbutylpiperidines such as fluspirilene, penfluridol, and pimozide, haloperidol decanoate, and indolones. Non-limiting examples of antianxiety-sedative agents include benzodiazepines such as chlordiazepoxide, diazepam, oxazepam, clorazepate, lorazepam, prazepam, alprazolam, and halazepam). Antidepressants or mood elevating agents include norepinephrine-reuptake inhibitors such as the tertiary amine tricyclics (e.g., amitriptyline, clomipramine, doxepin, and imipramine) and secondary amine tricyclics (e.g., amoxapine, desipramine, maprotiline, protriptyline, and nortriptyline); selective serotonin-reuptake inhibitors (SSRIs) such as fluoxetine, fluvoxamine, paroxetine, sertraline, citalopram, escitalopram, and venlafaxine; atypical antidepressants such as bupropion, nefazodone, and trazodone; noradrenergic and specific serotonergic antidepressants such as mirtazapine; and monoamine oxidase inhibitors such as phenelzine, tranylcypromine, and selegiline.

In certain ethnic groups as many as 10% of the adolescent population have a 2D6 haplotype that is associated with poor metabolism of many antidepressant medications. See Wong et al. (2001) Ann. Acad. Med. Singapore 29:401-406. Clinical genomic testing of these individuals has clear implications for their treatment and prognosis. In extreme cases, children who were poor metabolizers and who were not identified have had tragic outcomes. These negative case reports have included a reported death of a nine-year-old boy who was not recognized to be a poor 2D6 metabolizer. The treatment of this child with fluoxetine continued despite the development of multiple symptoms because these symptoms were not recognized as being related to his extremely high serum levels of fluoxetine. Sallee et al. (2000) J. Child Adol. Psychiatry 10(1):27-34. While careful clinical surveillance to identify unexpected side effects during the administration of low doses of medication is an alternative clinical strategy to genotyping, genomic testing of a plurality of genes encoding drug metabolizing enzymes (e.g. cytochrome P450 genes) and other target genes (e.g., genes involved in neurotransmission for psychotropic drugs) provides a safe method by which potentially dangerous side effects can be avoided in an affected patient.

Panels of Genes

The method includes obtaining a biological sample from a patient and obtaining the patient's genotype for a panel of genes. Typically, the panel of genes that are genotyped includes at least three cytochrome P450 genes. The cytochrome P450 genes can be selected from the P450 genes listed in Table 1. For example, the at least three cytochrome P450 genes can encode CYP2D6, 1A2, and 2C19. In embodiments for selecting antidepressant medications, the genotype of four cytochrome P450 genes (e.g., genes encoding CYP2D6, 2C19, 3A4, and 1A2) can be obtained. In other embodiments, the genotype of at least 5 cytochrome P450 genes (e.g., genes encoding CYP1A1, 1A2, 2D6, 2C19, and 3A4) can be obtained. In still other embodiments, the genotype of at least 10 cytochrome P450 genes (e.g., genes encoding CYP1A1, 1A2, 1B1, 2A6, 2B6, 208, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, and 3A5) can be obtained. Alleles for each of these cytochrome P450 genes are set forth in Table 1.

Substrates of CYP2D6 typically are weak bases with the cationic binding site located away from the carbon atom to be oxidized. In particular, substrates of CYP2D6 include amitriptyline, nortriptyline, haloperidol, and desipramine. Some individuals have altered CYP2D6 gene sequences that result in synthesis of enzymes devoid of catalytic activity or in enzymes with diminished catalytic activity. These individuals metabolize SSRIs and tricyclic antidepressants (TCAs) poorly. Duplication of the functional CYP2D6 gene also has been observed and results in ultrarapid metabolism of SSRIs and other drugs. Individuals without inactivating polymorphisms, deletions, or duplications have the phenotype of an extensive drug metabolizer and are designated as CYP2D6*1. The CYP2D6*2 allele has decreased enzymatic activity resulting from amino acid substitutions. The CYP2D6*3 and *4 alleles account for nearly 70% of the total deficiencies that result in the poor metabolizer phenotype. The polymorphism responsible for CYP2D6*3 (2549A>del) produces a frame-shift in the nRNA. A polymorphism involved with the CYP2D6*4 allele (1846G>A) disrupts mRNA splicing. These changes produce truncated forms of CYP2D6 devoid of catalytic activity. Other poor metabolizers are CYP2D6*5, *10, and *17. CYP2D6*5 is due to complete gene deletion. The polymorphisms in CYP2D6*10 and *17 produce amino acid substitutions in the CYP2D6 enzyme which have decreased enzyme activity. All of these polymorphisms are autosomal recessive. Consequently, only individuals who are homozygous or who are compound heterozygous for these polymorphisms are poor metabolizers. Individuals who are heterozygous, with one normal gene and one polymorphic gene, will have metabolism intermediate between the extensive (normal) and poor metabolizers.

CYP1A2 metabolizes many aromatic and heterocyclic amines including clozapine and imipramine. The CYP1A2*1F allele can result in a product with higher inducibility or increased activity. See Sachse et al. (1999) Br. J. Chin. Pharmacol. 47: 445-449. CYP2C19 also metabolizes many substrates including imipramine, citalopram, and diazepam. The CYP2C19*2A, *2B, *3, *4, *5A, *5B, *6, *7, and *8 alleles encode products with little or no activity. See Ibeanu et al. (1999) J. Pharmacol. Exp. Ther. 290: 635-640.

CYP1A1 can be associated with toxic or allergic reactions by extrahepatic generation of reactive metabolites. CYP3A4 metabolizes a variety of substrates including alprazolam. CYP1B1 can be associated with toxic or allergic reactions by extrahepatic generation of reactive metabolites and also metabolizes steroid hormones (e.g., 17β-estradiol). Substrates for CYP2A6 and CYP2B6 include valproic acid and bupropion, respectively. Substrates for CYP2C9 include Tylenol and antabuse (disulfuram). Substrates for CYP2E1 include phenytoin and carbamazepine. Decreases in activity in one or more of the cytochrome P450 enzymes can impact one or more of the other cytochrome P450 enzymes.

TABLE 1 Cytochrome P450 genes Cytochrome P450Gene Allele Polymorphism 1A1  *1A None *2 A2455G *3 T3205C *4 C2453A 1A2  *1A None  *1F −164C > A *3 G1042A 1B1 *1 None *2 R48G *3 L432V *4 N453S *11  V57C *14  E281X *18  G365W *19  P379L *20  E387K *25  R469W 2A6  *1A None   *1B CYP2A7 translocated to 3′-end *2 T479A *5 *1B + G6440T 2B6 *1 *2 R22C *3 S259C *4 K262R *5 R487C *6 Q172H; K262R *7 Q172H; K262R; R487C 2C8  *1A None   *1B −271C > A   *1C −370T > G *2 I269F *3 R139K; K399R *4 I264M 2C9 *1 None *2 R144C *3 I359L *5 D360E 2C18 m1 T204A m2 A460T 2C19  *1A None   *1B I331V  *2A Splicing defect   *2B Splicing defect; E92D *3 New stop codon 636G > A *4 GTG initiation codon, 1A > G *5(A, B) 1297C > T, amino acid change (R433W) *6 395G > A, amino acid change (R132Q) *7 IVS5 + 2T > A, splicing defect *8 358T > C, amino acid change (W120R) 2D6  *1A None *2 G1661C, C2850T   *2N Gene duplication *3 A2549 deletion *4 G1846A *5 Gene deletion *6 T1707 deletion *7 A2935C *8 G1758T *10  C100T *12  G124A *17  C1023T, C2850T *35  G31A 2E1  *1A None *1C, *1D (6 or 8 bp repeats) *2 G1132A *4 G476A *5 G(−1293)C *5 C(−1053)T *7 T(−333)A *7 G(−71)T *7 A(−353)G 3A4  *1A None   *1B A(−392)G *2 Amino acid change (S222P) *5 Amino acid change (P218R) *6 Frameshift, 831 ins A *12  Amino acid change (L373F) *13  Amino acid change (P416L)  *15A Amino acid change (R162Q) *17  Amino acid change (F189S, decreased)  *18A Amino acid change (L293P, increased) 3A5  *1A None *3 A6986G *5 T12952C *6 G14960A

Typically, to select a medication, the genotype of other target genes also is obtained in addition to the genotype of the genes encoding drug-metabolizing enzymes. Target genes can encode products that relate to the ability of the patient to respond to a particular class of medication. For example, to select an antidepressant, the target genes can be the serotonin transporter gene and serotonin receptor 2A gene. To select an antipsychotic, the target gene can be a dopamine transporter gene.

Table 2 sets forth the alleles for target genes that can be genotyped in addition to the genes encoding drug-metabolizing enzymes. For example, the panel of genes to be genotyped can include one or more dopamine receptor genes (e.g., genes encoding dopamine receptors D1, D2, D3, D4, D5, and/or D6), a dopamine transporter gene, one or more serotonin receptor genes (e.g., genes encoding serotonin receptors 1A, 1B, 1D, 2A, or 2C), a catechol-O-methyl transferase (COMT) gene, or a tryptophan hydroxylase gene. In one embodiment, a COMT gene and a tryptophan hydroxylase gene are included on the panel with cytochrome P450 genes. In other embodiments, one or more dopamine receptor genes, one or more serotonin receptor genes, a COMT gene, and a tryptophan hydroxylase gene are assessed in combination with the cytochrome P450 genes.

TABLE 2 Gene Symbol Polymorphism Dopamine Transporter DAT1, 40 bp VNTR SLC6A3 10 repeat allele G710A, Q237R C124T, L42F Dopamine Receptor D1 DRD1 DRD1 B2 T244G C179T G127A T11G C81T T595G, S199A G150T, R50S C110G, T37R A109C, T37P Dopamine Receptor D2 DRD2 TaqI A A1051G, T35A C932G, S311C C928, P310S G460A, V154I Dopamine Receptor D3 DRD3 BalI in exon I MspI DRD3 1 Gly/Ser (allele 2) A25G, S9G Dopamine Receptor D4 DRD4 48 repeat in exon 3 7 repeat allele 12/13 bp insertion/deletion T581G, V194G C841G, P281A Dopamine Receptor D5 DRD5 T978C L88F A889C, T297P G1252A, V418I G181A, V61M G185C, C62S T263G, R88L G1354A, W455 Tryptophan TPH A218C Hydroxylase A779C G-5806T A-6526G (CT)_(m)(CA)_(n)(CT)_(p) allele 194 in 3′ UTR, 5657 bp distant from exon 11 Serotonin Transporter 5-HTTR Promoter repeat (44 bp insertion (L)/deletion(S) (L = Long form; S = Short form) Exon 2 variable repeat A1815C G603C G167C Serotonin Receptor 1A HTR1A RsaI G815A, G272D G656T, R219L C548T, P551L A82G, I28V G64A, G22S C47T, P16L Serotonin Receptor 1B HTR1B G861C G861C, V287V T371G, F124C T655C, F219L A1099G, I367V G1120A, E374K Serotonin Receptor 1D HTR1D G506T C173T C794T, S265L Serotonin Receptor 2A HTR2A C74A T102C T516C C1340T C1354T Serotonin Receptor 2C HTR2C G796C C10G, L4V G68C, C23S Catechol-o- COMT G158A (Also known methyltransferase as Val/Met) G214T A72S G101C C34S G473A

For example, to select an antidepressant, the genotype of four cytochrome P450 genes (e.g., genes encoding 2D6, 2C19, 3A4, and 1A2), the serotonin transporter (5HTTR) gene, and serotonin receptor 2A (HTR2A) gene can be assessed in a patient. In particular, it can be determined if the patient contains the CYP2D6*2, *3, *4, *10, *17, or *5 del allele, the CYP2C19*2(A,B), *3, *4, *5 (A,B), *6, *7, *8 alleles, the CYP3A4*1B, *2, *5, *6, *12, *13, *15A, *17, or *18A alleles, the CYP1A2*1F allele, the short or long form of the serotonin transporter gene, and the HTR2A T102C polymorphism. Each of these genes influences the metabolism of at least one antidepressant medication or is associated with a better treatment response. As described herein, an algorithm has been created based on a set of six rules relating to the genotype of these six genes. Based on this algorithm, medication profiles are provided for a given patient based on the patient's genotype, allowing a clinician to select an acceptable antidepressant without the trial and error of determining if the patient will respond or tolerate a particular antidepressant.

Determining Genotype

Genomic DNA generally is used to determine genotype, although mRNA also can be used. Genomic DNA is typically extracted from a biological sample such as a peripheral blood sample, but can be extracted from other biological samples, including tissues (e.g., mucosal scrapings of the lining of the mouth or from renal or hepatic tissue). Routine methods can be used to extract genomic DNA from a blood or tissue sample, including, for example, phenol extraction. Alternatively, genomic DNA can be extracted with kits such as the QIAamp® Tissue Kit (Qiagen, Chatsworth, Calif.), Wizard® Genomic DNA purification kit (Promega) and the A.S.A.P.™ Genomic DNA isolation kit (Boehringer Mannheim, Indianapolis, Ind.).

Typically, an amplification step is performed before proceeding with the genotyping. For example, polymerase chain reaction (PCR) techniques can be used to obtain amplification products from the patient. PCR refers to a procedure or technique in which target nucleic acids are enzymatically amplified. Sequence information from the ends of the region of interest or beyond typically is employed to design oligonucleotide primers that are identical in sequence to opposite strands of the template to be amplified. PCR can be used to amplify specific sequences from DNA as well as RNA, including sequences from total genomic DNA or total cellular RNA. Primers are typically 14 to 40 nucleotides in length, but can range from 10 nucleotides to hundreds of nucleotides in length. General PCR techniques are described, for example in PCR Primer: A Laboratory Manual, Ed, by Dieffenbach, C. and Dveksler, G., Cold Spring Harbor Laboratory Press, 1995. When using RNA as a source of template, reverse transcriptase can be used to synthesize complementary DNA (cDNA) strands. Ligase chain reaction, strand displacement amplification, self-sustained sequence replication or nucleic acid sequence-based amplification also can be used to obtain isolated nucleic acids. See, for example, Lewis (1992) Genetic Engineering News 12(9):1; Guatelli et al. (1990) Proc. Natl. Acad. Sci. USA 87:1874-1878; and Weiss (1991) Science 254:1292-1293.

Primers typically are single-stranded or double-stranded oligonucleotides that are 10 to 50 nucleotides in length, and when combined with mammalian genomic DNA and subjected to PCR conditions, is capable of being extended to produce a nucleic acid product corresponding to a region of interest within a gene. Typically, PCR products are at least 30 nucleotides in length (e.g., 30, 35, 50, 100, 250, 500, 1000, 1500, or 2000 or more nucleotides in length). Primers such as those listed in Table 7 are particularly useful for producing PCR products for the genes encoding the dopamine transporter, dopamine receptors, tryptophan hydroxylase, serotonin transporter, serotonin receptors, and COMT. Specific regions of mammalian DNA can be amplified (i.e., replicated such that multiple exact copies are produced) when a pair of oligonucleotide primers is used in the same PCR reaction, wherein one primer contains a nucleotide sequence from the coding strand of a nucleic acid and the other primer contains a nucleotide sequence from the non-coding strand of the nucleic acid. The “coding strand” of a nucleic acid is the nontranscribed strand, which has the same nucleotide sequence as the specified RNA transcript (with the exception that the RNA transcript contains uracil in place of thymidine residues), while the “non-coding strand” of a nucleic acid is the strand that serves as the template for transcription.

A single PCR reaction mixture may contain one pair of oligonucleotide primers. Alternatively, a single reaction mixture may contain a plurality of oligonucleotide primer pairs, in which case multiple PCR products can be generated (e.g., 5, 10, 15, or 20 primer pairs). Each primer pair can amplify, for example, one exon or a portion of one exon. Intron sequences also can be amplified.

Exons or introns of a gene of interest can be amplified then directly sequenced. Dye primer sequencing can be used to increase the accuracy of detecting heterozygous samples. Alternatively, one or more of the techniques described below can be used to determine genotype.

For example, allele specific hybridization can be used to detect sequence variants, including complete haplotypes of a mammal. See, Stoneking et al., 1991, Am. J. Hum. Genet. 48:370-382; and Prince et al., 2001, Genome Res., 11(1):152-162. In practice, samples of DNA or RNA from one or more mammals can be amplified using pairs of primers and the resulting amplification products can be immobilized on a substrate (e.g., in discrete regions). Hybridization conditions are selected such that a nucleic acid probe can specifically bind to the sequence of interest, e.g., the variant nucleic acid sequence. Such hybridizations typically are performed under high stringency as some sequence variants include only a single nucleotide difference. High stringency conditions can include the use of low ionic strength solutions and high temperatures for washing. For example, nucleic acid molecules can be hybridized at 42° C. in 2×SSC (0.3M NaCl/0.03 M sodium citrate/0.1% sodium dodecyl sulfate (SDS) and washed in 0.1×SSC (0.015M NaCl/0.0015 M sodium citrate), 0.1% SDS at 65° C. Hybridization conditions can be adjusted to account for unique features of the nucleic acid molecule, including length and sequence composition. Probes can be labeled (e.g., fluorescently) to facilitate detection. In some embodiments, one of the primers used in the amplification reaction is biotinylated (e.g., 5′ end of reverse primer) and the resulting biotinylated amplification product is immobilized on an avidin or streptavidin coated substrate (e.g., in discrete regions).

Allele-specific restriction digests can be performed in the following manner. For nucleotide sequence variants that introduce a restriction site, restriction digest with the particular restriction enzyme can differentiate the alleles. For sequence variants that do not alter a common restriction site, mutagenic primers can be designed that introduce a restriction site when the variant allele is present or when the wild type allele is present. A portion of the nucleic acid of interest can be amplified using the mutagenic primer and a wild type primer, followed by digest with the appropriate restriction endonuclease.

Certain variants, such as insertions or deletions of one or more nucleotides, change the size of the DNA fragment encompassing the variant. The insertion or deletion of nucleotides can be assessed by amplifying the region encompassing the variant and determining the size of the amplified products in comparison with size standards. For example, a region of a gene of interest can be amplified using a primer set from either side of the variant. One of the primers is typically labeled, for example, with a fluorescent moiety, to facilitate sizing. The amplified products can be electrophoresed through acrylamide gels with a set of size standards that are labeled with a fluorescent moiety that differs from the primer.

PCR conditions and primers can be developed that amplify a product only when the variant allele is present or only when the wild type allele is present (MSPCR or allele-specific PCR). For example, patient DNA and a control can be amplified separately using either a wild type primer or a primer specific for the variant allele. Each set of reactions is then examined for the presence of amplification products using standard methods to visualize the DNA. For example, the reactions can be electrophoresed through an agarose gel and the DNA visualized by staining with ethidium bromide or other DNA intercalating dye. In DNA samples from heterozygous patients, reaction products would be detected in each reaction. Patient samples containing solely the wild type allele would have amplification products only in the reaction using the wild type primer. Similarly, patient samples containing solely the variant allele would have amplification products only in the reaction using the variant primer, Allele-specific PCR also can be performed using allele-specific primers that introduce priming sites for two universal energy-transfer-labeled primers (e.g., one primer labeled with a green dye such as fluoroscein and one primer labeled with a red dye such as sulforhodarnine). Amplification products can be analyzed for green and red fluorescence in a plate reader. See, Myakishev et al., 2001, Genome 11(1):163-169.

Mismatch cleavage methods also can be used to detect differing sequences by PCR amplification, followed by hybridization with the wild type sequence and cleavage at points of mismatch. Chemical reagents, such as carbodiimide or hydroxylamine and osmium tetroxide can be used to modify mismatched nucleotides to facilitate cleavage.

Kits also are available commercially to detect many of the cytochrome P450 variants. For example, TAG-IT™ kits are available from Tm Biosciences Corporation (Toronto, Ontario).

Selecting Medications

After the genotype is determined for each gene on the panel, the medication can be selected. Typically, selecting includes correlating genotype of the cytochrome P450 genes with capacity of each cytochrome P450 enzyme encoded by each cytochrome P450 gene to metabolize the medication. The genotype of other target genes on the panel, e.g., the serotonin transporter and/or dopamine transporter genes, can be correlated with the ability of the patient to respond to the medication.

An algorithm can be used to select the most appropriate medications for an individual patient. The design of the algorithm requires the initial identification of the phenotype, which provides a preliminary identification of the universe of possible medications. At the next step of the algorithm, the results of the target gene analyses can be sequentially entered. The potential list of appropriate medications can be subsequently rank-ordered based on specific co-factors in the algorithmic equation, which assign a positive, negative, or neutral probability score to each of the identified medications in the originally identified set of possible choices. This process adjusts the rank ordering based on the genotypic polymorphism carried by the patient. At the next entry point in the algorithmic equation, the results of the CYP gene analyses can be introduced. This algorithmic analysis is designed to place the medications in three categories: 1) medications that are acceptable for use, i.e., the medication has a high probability of normal metabolism within an individual having a particular genotype, 2) medications that can be used with caution (e.g., medication may require some dosing adjustment based on atypical metabolism); and 3) medications that should be avoided or used with caution and monitoring, e.g., due to potential difficulties in dosing. At this point in the selection process, data related to the medication response of first- and second-degree relatives of the patient can be entered into the algorithmic equation, which pertains to the medication selection of drugs in the first category that has been identified. An adjustment of the rank-ordered, appropriate medications then can be calculated based on clinical responses by family members.

The selection of an appropriate medication further can be enhanced by including both target data and data related to drug metabolism. This can determine the impact of the CYP products on the clinical response of a particular patient. For example, inclusion of target data and data related to drug metabolism provides the amount of available drug, the ability of the patient to utilize the drug, and information about the quality of the receptor target of the drug, providing a rational approach to selection of medication.

An example of this process would be the selection of an appropriate antidepressant for a given patient. Once the phenotype of depression is established, the initial universe of antidepressant medications is identified. For example, citalopram, fluvoxamine, bupropion, escitalopram, sertraline, mirtazapine, flouxetine, venlafaxine, amitriptyline, imipramine, paroxetine, and nortriptyline can be the initial universe of antidepressant medications. It is understood that this initial universe of medications can change as additional antidepressants become available or as additional pharmacogenomic data are available. Subsequently, the data from the target genotyping is entered. If during the first components of the algorithm, it is revealed that the genotype reflects poor metabolism for 2D6 and normal metabolism for 1A2, the rank ordering of medication would be to identify fluvoxamine and buproprion as medications that are likely to be both effective and appropriately managed at recommended doses. At the next stage of the algorithm, if an individual is homozygous for the long form of the serotonin transporter gene, fluvoxamine would be confirmed as being acceptable for use. Output from the algorithm then can be integrated with historical data. For example, if a family member had responded well to fluvoxamine, this would confirm that the medication is acceptable for use, or, if a first or second degree relative had a problematic response to this medication, an alternative could be chosen.

Computer Systems

Techniques described herein can be implemented in a computer system having a processor that executes specific instructions in a computer program. The computer system may be arranged to output a medication profile based on receiving a patient's genotype. Particularly, the computer program may include instructions for the system to select the most appropriate medication (e.g., a psychotropic medication) for an individual patient.

The following are examples of features that may be included in a system. The computer program may be configured such that the computer system can identify the phenotype based on received data and provide a preliminary identification of the universe of possible medications. The system may be able to rank-order the identified medications based on specific co-factors in the algorithmic equation. The system may be able to adjust the rank ordering based on the genotypic polymorphism(s) carried by the patient. The system may be able to adjust the rank ordering based on clinical responses, such as by family members of the patient.

FIG. 1 is a block diagram of a computer system 100 that can be used in the operations described above, according to one embodiment. The system 100 includes a processor 110, a memory 120, a storage device 130 and an input/output device 140. Each of the components 110, 120, 130 and 140 are interconnected using a system bus 150. The system may include analyzing equipment 160 for determining the patient's genotype.

The processor 110 is capable of processing instructions for execution within the system 100. In one embodiment, the processor 110 is a single-threaded processor. In another embodiment, the processor 110 is a multi-threaded processor. The processor 110 is capable of processing instructions stored in the memory 120 or on the storage device 130, including for receiving or sending information through the input/output device 140.

The memory 120 stores information within the system 100. In one embodiment, the memory 120 is a computer-readable medium. In one embodiment, the memory 120 is a volatile memory unit. In another embodiment, the memory 120 is a non-volatile memory unit.

The storage device 130 is capable of providing mass storage for the system 100. In one embodiment, the storage device 130 is a computer-readable medium. In various different embodiments, the storage device 130 may be a floppy disk device, a hard disk device, an optical disk device, or a tape device.

The input/output device 140 provides input/output operations for the system 100. In one embodiment, the input/output device 140 includes a keyboard and/or pointing device. In one embodiment, the input/output device 140 includes a display unit for displaying graphical user interfaces.

The system 100 can be used to build a database. FIG. 2 shows a flow chart of a method 200 for building a database for use in selecting a medication for a patient. Preferably, the method 200 is performed in the system 100. For example, a computer program product can include instructions that cause the processor 110 to perform the steps of the method 200. The method 200 includes the following steps.

Receiving, in step 210, a plurality of genotypes 170 for a panel of genes. A computer program in the system 100 may include instructions for presenting a suitable graphical user interface on input/output device 140, and the graphical user interface may prompt the user to enter the genotypes 170 using the input/output device 140, such as a keyboard.

Receiving, in step 220, a plurality of medication profiles 180. The medication profiles 180 are specified based on the genotypes 170. The user may enter the medication profiles 180 using the input/output device 140, such as a keyboard. For example, the medication profile 180 may include information 190 regarding at least one medication.

Storing, in step 230, the received genotypes 170 and the medication profiles 180 such that each medication profile 180 is associated with one of the genotypes 170. The system 100 may store the medication profiles 180 and the genotypes 170 in the storage device 130. For example, when the storing is complete, the system 100 can identify a particular one of the medication profiles 180 that is associated with a specific genotype 170. Having identified the medication profile 180, the system 100 can access the information 190 contained within the identified medication profile 180, as will be described in the following example.

The system 100 may be used for selecting a medication. FIG. 3 shows a flow chart of a method 300 of selecting a medication for a patient. Preferably, the method 300 is performed in the system 100. For example, a computer program product can include instructions that cause the processor 110 to perform the steps of the method 300. The method 300 includes the following steps.

Receiving, in step 310, a patient's genotype for a panel of genes. The genotype may be entered by a user via input/output device 140. For example, the user may obtain the patient's genotype for a panel of genes using the analyzing equipment 160 (which may or may not be connected to the system 100). The user may type the patient's genotype on input/output device 140, such as a keyboard, for receipt by the system 100.

The genotype may be received directly from the analyzing equipment 160. For example, analyzing equipment 160 may include a processor and suitable software such that it can communicate over a network. The system 100 may be connected to the analyzing equipment 160 through input/output device 140, such as a network adapter, and directly receive the patient's genotype.

Identifying, in step 320, one of the medication profiles 180 that is associated with the patient's genotype. For example, the system 100 may perform a database search in the storage device 130. Particularly, the system 100 may access the genotype 170 for individual medication profiles 180 until a match is found. Optional step 325 will be described below.

Outputting, in step 330, the identified medication profile 180 in response to receiving the patient's genotype. The system may output the identified medication profile 180 through input/output device 140. For example, the identified medication profile may be printed or displayed in a suitable graphical user interface on a display device. As another example, the system 100 may transmit the identified medication profile over a network, such as a local area network or the Internet, to which the input/output device 140 is connected.

The medication profiles 180 can be created such that there is flexibility in how the system 100 outputs them. For example, the information 190 in one or more of the medication profiles 180 may include a ranking of several medications. The program may include instructions for applying rules to the received patient's genotype and adjust the ranking accordingly. In such implementations, the method 300 may include optional step 325 of adjusting the ranking before outputting the identified medication profile. For example, the system 100 may receive a genotypic polymorphism carried by the patient (optionally in the same way the patient's genotype was received) and adjust the ranking accordingly in step 325. As another example, step 325 may involve adjusting the ranking based on a clinical response. The clinical response may be received by the system 100 in the same way as the patient's genotype. For example, the ranking can be adjusted based on a clinical response by a member of the patient's family.

The medication profiles 180 may be updated as necessary. For example, the introduction of a new medication on the market may prompt a revision of one or more existing medication profiles. A new medication may also be the basis for creating a new medication profile. The adjustment or creation of medication profiles may be done substantially as described above.

The medication profiles 180 may be used for medication selection in the same system where they were created; or in a different system. That is, the system 100 may first be used for building a database of the medication profiles 180, and the system 100 may thereafter be used to select a medication profile for the genotype of a specific patient. As another example, one or more medication profiles 180 may be transmitted within a computer readable medium such as a global computer network for remote processing according to the invention.

Articles of Manufacture

In one embodiment, the article of manufacture is a composition containing oligonucleotide primers. Typically, such a composition will contain a first oligonucleotide primer and a second oligonucleotide primer, each 10 to 50 nucleotides in length, which can be combined with genomic DNA from a mammal and subjected to PCR conditions, to produce a nucleic acid product that corresponds to a region of interest within a target gene. A composition also may contain buffers and other reagents necessary for PCR (e.g., DNA polymerase or nucleotides). Furthermore, a composition may contain one or more additional pairs of oligonucleotide primers (e.g., 5, 10, 15, or 20 primer pairs), such that multiple nucleic acid products can be generated.

In other embodiments, articles of manufacture include populations of nucleic acid molecules immobilized on a substrate. Suitable substrates provide a base for the immobilization of the nucleic acids, and in some embodiments, allow immobilization of nucleic acids into discrete regions. In embodiments in which the substrate includes a plurality of discrete regions, different populations of isolated nucleic acids can be immobilized in each discrete region. The different populations of nucleic acid molecules independently can include nucleic acid molecules for detecting one or more of the alleles set forth in Tables 1 and 2.

Suitable substrates can be of any shape or form and can be constructed from, for example, glass, silicon, metal, plastic, cellulose or a composite. For example, a suitable substrate can include a multiwell plate or membrane, a glass slide, a chip, or polystyrene or magnetic beads. Nucleic acid molecules or polypeptides can be synthesized in situ, immobilized directly on the substrate, or immobilized aria a linker, including by covalent, ionic, or physical linkage. Linkers for immobilizing nucleic acids and polypeptides, including reversible or cleavable linkers, are known in the art. See, for example, U.S. Pat. No. 5,451,683 and WO98/20019. Immobilized nucleic acid molecules typically are about 20 nucleotides in length, but can vary from about 10 nucleotides to about 1000 or more nucleotides in length.

The invention will be further described in the following examples, which do not limit the scope of the invention described in the claims.

EXAMPLES Example 1 CYP450 Genotype of Patients

This experiment highlights the results of an assessment of adolescents with atypical responses to selective serotonin reuptake inhibitors. The following CYP2D6 alleles that produce deficiencies were detected: *2 (1661G>C and 2850C>T), *3 (2549A>del), *4 (16610>C and 1846G>A), *10 (100C>T, 16610>C and 1846G), *17 (1023C>T and 2850C>T). Deletion of CYP2D6 (*5 del) was detected by electrophoresis following PCR amplification. Tandem duplication of the CYP2D6 gene was detected by amplification of a duplication-specific sequence with detection following electrophoresis.

Genomic DNA was amplified by the polymerase chain reaction in an Applied Biosystems Thermal Cycler using the primers listed below. The lower case letters “f” and “r” imply the direction of synthesis (“f” indicates a forward/sense synthesis reaction and “r” indicates a reverse/antisense synthesis reaction). The amplified DNA was hybridized with oligonucleotide probes (described below) specific for the normal and polymorphic variants and detected using the Nanogen Molecular Biology Workstation. The oligonucleotide probes were fluorescently labeled (e.g., CY3 or CY5).

PreNest Oligonucleotides (Primary Amplification)

Amplification Primers:

CYP2D6 Prenest B F (Bases 1281-1302, SEQ ID NO: 1) 5′-CCA GAA GGC TTT GCA GGC TTC A-3′ CYP2D6 Prenest B R (Bases 6357-6378, SEQ ID NO: 2) 5′-ACT GAG CCC TGG GAG GTA GGT A-3′

Nest Oligonucleotides (Secondary Amplification)

CYP2D6 (C100T)

Amplification Primers:

CYP2D6 (100)B f (Bases 1579-1594, SEQ ID NO: 3) 5′GGC CTA CCC TGG GTA AGG GCC TGG AGC AGG A-3′ CYP2D6 (100)B r (Bases 2394-2411, SEQ ID NO: 4) 5′-/5Bio/CCT GGT CGA AGC AGT AT-3′

Probes and Stabilizers

CYP2D6 C100T wild type (Bases 1700-1710, SEQ ID NO: 5) 5′-/5Cy3/GCA CGC TAC C-3′ CYP2D6 C100T polymorphism (Bases 1700-1710, SEQ ID NO: 6) 5′-/5Cy5/GCA CGC TAC T-3′ CYP2D6 C100T stabilizer (Bases 1711-1731, SEQ ID NO: 7) 5′CAC CAG GCC CCC TGC CAC TG-3′

CYP2D6 (C1023T)

Amplification Primers:

CYP2D6 (1023)B f (Bases 2394-2413, SEQ ID NO: 8) 5′-CCT GCT CAC TCC TGG TAG CC-3′ CYP2D6(1023) B r (Bases 2760-2779, SEQ ID NO: 9) 5′-/5Bio/CTG TTT CAT GTC CAC GAC C-3′

Probes and Stabilizers

CYP2D6 (1023) wild type (Bases 2632-2642, SEQ ID NO: 10) 5′-/5Cy3/TGC CCA TCA C-3′ CYP2D6 (1023) polymorphism (Bases 2632-2642, SEQ ID NO: 11) 5′-/5Cy3/TGC CCA TCA T-3′ CYP2D6 (1023) stabilizer (Bases 2643-2667, SEQ ID NO: 12) 5′CCA GAT CCT GGG TTT CGG GCC GCG-3′

CYP2D6 (01661C and G1846A)

Note: There are 2 polymorphisms on this amplicon.

Amplification Primers

CYP2D6 (1661, 1846)B f (Bases 3266-3282, SEQ ID NO: 13) 5′-CAG AGG CGC TTC TCC G-3′ CYP2D6 (1661, 1846)B r (Bases 3630-3650, SEQ ID NO: 14) 5′-/5Bio/CTC GGT CTC TCG CTC CGC AC-3′

Probes and Stabilizers

CYP2D6(1661) wild type (Bases 3266-3276, SEQ ID NO: 15) 5′-/5Cy3/CTT CTC CGT G-3′ CYP2D6(1661) polymorphism (Bases 3266-3276, SEQ ID NO: 16) 5′-/5Cy3/CTT CTC CGT C-3′ CYP2D6(1661) stabilizer (Bases 3277-3303, SEQ ID NO: 17) 5′-TCC ACC TTG CGC AAC TTG GGC CTG GG-3′ CYP2D6(1846) wild type (Bases 3460-3470, SEQ ID NO: 18) 5′-/5Cy3/CCA CCC CCA G-3′ CYP2D6(1846) polymorphism (Bases 3460-3470, SEQ ID NO: 19) 5′-/5Cy3/CCA CCC CCA A-3′ CYP2D6(1846) stabilizer (Bases 3471-3492, SEQ ID NO: 20) 5′-GAC GCC CCT TTC GCC CCA ACG-3′

CYP2D6 (A2549del and C2850T)

Amplification Primers:

CYP2D6 (A2549del and C2850T)B f (Bases 4001-4022, SEQ ID NO: 21) 5′-TGA GAC TTG TCC AGG TGA AC-3′ CYP2D6(A2549del and C2850T)B r (Bases 4558-4578, SEQ ID NO: 22) 5′-/5Bio/CCC AGA TGG GCT CAC GCT GC-3′

Probes and Stabilizers

CYP2D6 (2549) wild type (Bases 4161-4171, SEQ ID NO: 23) 5′-/5Cy3/ACT GAG CAC A-3′ CYP2D6 (2549) polymorphism (Bases 4161-4171, SEQ ID NO: 24) 5′-/5Cy3/ACT GAG CAC G-3′ CYP2D6 (2549) stabilizer (Bases 4172-4199, SEQ ID NO: 25) 5′-GGA TGA CCT GGG ACC CAG CCC AGC C3′ CYP2D6 (2850) wild type (Bases 4462-4472, SEQ ID NO: 26) 5′-/5Cy3/GAG AAC CTG C-3′ CYP2D6 (2850) polymorphism (Bases 4462-4472, SEQ ID NO: 27) 5′-/5Cy3/GAG AAC CTG T-3′ CYP2D6 (2850) stabilizer (Bases 4473-4500, SEQ ID NO: 28) 5′-GCA TAG TGG TGG CTG ACC TGT TCT CTG-3′

PCR Amplification.

A. PreNest (Initial Amplification);

CYP2D6 PreNest B

The initial PreNest PCR for the CYP2D6(C1661G), CYP2D6 (G1846A), CYP2D6(A2549Del), CYP2D6(C100T), CYP2D6(C2850T), and, CYP2D6 (C1023T) polymorphisms was performed with the following reagents per sample: 0.6 μL 25 μM PreNest B F, 0.6 μL 25 μM PreNestB R, 2.5 μL 10 mM dNTP's, 39.55 μL dH₂O, 5.0 μL Expand Buffer 3 (Roche), 0.75 μL Expand Enzyme Mix (Roche), and 1.0 μl genomic DNA in a total volume of 50 μL. PCR was performed by holding at 94° C. for 2 minutes then 30 cycles of denaturation at 94° C. for 30 seconds, annealing at 67° C. for 30 seconds, and extension at 68° C. for 7 minutes, followed by a final extension at 72° C. for 10 minutes. The expected PCR product was ≈4.8-5 Kb.

All subsequent PCR's were performed using the PreNest B product as the template for the amplification.

CYP2D6(100) B

The (C100T) nested PCR was performed using 1 μl of the initial PreNest B amplicon as template. The product of this PCR was used to detect the CYP2D6(C100T) mutation and uses the following reagents (25 μL total volume): 1.0 μl 25 μM CYP2D6(100) B f, 1.0 μl 25 μM CYP2D6(100) B r, 12.5 μl Qiagen Hot Start Master Mix, and 9.5 μl dH₂O. PCR was performed by holding at 95° C. for 15 minutes, then 30 cycles of denaturation at 95° C. for 30 sec, annealing at 57° C. for 30 sec, and extension at 72° C. for 1 min, followed by a final extension at 72° C. for 7 min.

CYP2D6(1023) B

A (C1023T) nested PCR was performed using 1 μl of PreNest B product as template. The product of this PCR was used to detect the CYP2D6 (C1023T) mutation and uses the following reagents (25 μL total volume); 1.0 μl 25 μM CYP2D6 (1023) B f, 1.0 μl 25 μM CYP2D6 (1023)B r, 12.5 μl Amplitaq Gold Master Mix, and 9.5 μl dH₂O. PCR was performed by holding at 95° C. for 10 min, then 30 cycles of denaturation at 95° C. for 30 sec, annealing at 59° C. for 30 sec, and extension at 72° C. for 1 min, followed by a final extension at 72° C. for 7 min.

CYP2D6(C1661G and G1846A) B

The (C1661 G and G1846A) nested PCR was performed using 1 μl of PreNest B product as template. The product of this PCR was used to detect the CYP2D6 (C1661G) and CYP2D6 (G1846A) mutations and uses the following reagents (25 μL total volume): 1.0 μl 25 μM CYP2D6 (1661,1846)B f, 1.0 μl 25 μM CYP2D6 (1661,1846)B r, 12.5 μl 2× Amplitaq Gold Master Mix, and 9.5 μl dH₂O. PCR was performed by holding at 95° C. for 10 min, then 30 cycles of denaturation at 95° C. for 30 sec, annealing at 63° C. for 30 sec, and extension at 72° C. for 1 min, followed by a final extension at 72° C. for 7 min.

CYP2D6 (A2549Del and CYP2D6C2850T) B

The (A2549Del and C2850T) nested PCR was performed using 1 μl of PreNest B product as template. The product of this PCR was used to detect the CYP2D6 (A2549Del) and CYP2D6 (C2850T) mutations and uses the following reagents (25 μL total volume): 1.0 μl 25 μM CYP2D6 (2549,2850)B f, 1.0 μl 25 μM CYP2D6 (2549,2850)B r, 12.5 μl 2× Amplitaq Gold Master Mix, and 9.5 μl dH₂O. PCR was performed by holding at 95° C. for 10 min, then 30 cycles of denaturation at 95° C. for 30 sec, annealing at 57° C. for 30 sec, and extension at 72° C. for 1 min, followed by a final extension at 72° C. for 10 min.

Polymorphisms were detected using the Nanogen Molecular Biology Workstation according to the Nanogens User's manual. Samples were desalted on Millipore desalting plates. In a 1.5 ml centrifuge tube, the following reagents were added: 2.0 μl of wild type probe, 2.0 μl of polymorphism probe, 2.0 μl of stabilizer, and 44 μl high salt buffer. Each assay had probes and stabilizers that were labeled and specific to it. The reagents were vortexed, spun briefly, and dispensed over the Nanochip array for 5 minutes and placed in a dark container. The chips were read with a trending method from 24° C. to 50° C., with readings every two degrees. It was confirmed that the fluorescent probe hybridized at 24° C. and that all probes had de-hybridized at 50° C. If the relative fluorescence units (RFUs) were above 80 on any chip at 50° C., the chips were stripped in 200 μl of 1.0 M NaOH for 10 minutes and washed three times with dH₂O and once with high salt buffer. A heterozygote control (known heterozygote, one to “normalize” or scale the RFUs and the other to confirm the correcting trending pattern observed in previously genotyped heterozygotes) and a negative control (no DNA) were included for each polymorphism. For background subtraction, L-histidine was used.

Tables 3-5 summarize the total number of participants for different age groups and the CYP2D6 inferred metabolic phenotype detected using the above procedure. Table 6 provides the particular 2D6 haplotypes that were observed in the patients that had antidepressant medication resistant subtypes of depression.

TABLE 3 Total Number of Participants in Each Grouping (All Ages; N = 86) CYP2D6 Inferred Metabolic Phenotype Ultra- Poor Intermediate Extensive Rapid Response to Non-Response 5 16 6 0 Psychotropic Side-Effect 5 16 7 0 Both 8 17 5 1 Total 18 49 18 1

TABLE 4 Total Number of Participants in Each Grouping (Ages: <19; N = 12) CYP2D6 Inferred Metabolic Phenotype Ultra- Poor Intermediate Extensive Rapid Response Non-Response 0 3 0 0 to Side-Effect 3 3 0 0 Psychotropic Both 1 1 0 1 Total 4 7 0 1

TABLE 5 Total Number of Participants in Each Grouping (Ages: >18; N = 74) CYP2D6 Inferred Metabolic Phenotype Ultra- Poor Intermediate Extensive Rapid Response to Non-Response 5 13 6 0 Psychotropic Side-Effect 2 13 7 0 Medication Both 7 16 5 0 Total 14 42 18 0

TABLE 6 Antidepressant Medication Resistant Subtypes of Depression¹ 2D6 Enz Sub# 2D6 Haplo Act 100MR *1/*2P Hz*2 I 101MR Hz *4 I 102MR *1/*10 H2*10 I 103MR Hz *4 I 104MR *2/*2P Homo *2 I-P 105MR Hz *11 I 106MR Hz *4 I 107MR Homo *3 None 108MR Hz *4 I 109MR *2/*2 Homo *2 I-P 110MR Hz *2 I 111MR Hz *2 I 112MR WT N 113MR Hz *2 I 114MR WT N 115MR WT N 116MR *3/*4 Compd Hz *3/*4 P 117MR Hz *4 I 118MR Hz *2 I 119MR Compound Hz *2/*4 P 120MR Hz *4 I 121MR Hz *4 I 122MR Hz *4 I 123MR WT N 124MR Hz *12 I 125MR Hz *4 I 126MR Hz *4 I 128MR WT N 129MR Hz *4 I 130MR Hz *1E N 131MR Homo *4 P 132MR Hz *4 I 133MR Hz *2 I 134MR Compound Hz *2/*3 P 135MR Compound Hz *3/*4 P 136MR Hz * I 137MR Hz *2 I 138MR Compound Hz *2/*4 P 139MR Compound Hz *2/*5 P 141MR Homo *5 None 142MR Compound Hz *4/*15 P 143MR Hz *4 I 144MR Hz *17 or *34 I 145MR Hz *3 I 146MR Compound Hz *4/*5 P 147MR WT N 148MR Hz *2 I 149MR WT N 150MR Hz *2 I 155MR WT N 156MR Hz *4 I 157MR Hz *2 I 158MR Hz *4 I 159MR Hz *2 I 160MR WT N 161MR Hz *2 I 162MR WT N 164MR Hz *2 I 165MR Hz *2 I 166MR WT N 167MR Compound Hz *2/*4 P 168MR Compound Hz *2/*4 P 169MR WT N 170MR Hz *2 I 171MR Compound Hz *2/*3 P 172MR Hz *2 w/ duplication Inc. act. 173MR Hz *2 I 174MR Hz *4 w/ balanced duplication N 175MR Hz *2 I 176MR Hz *4 I 177MR Hz *2 I 179MR Homo *2 I-P 180MR Homo *2 I-P 181MR WT N 182MR Hz *4 I 183MR WT N 184MR Hz *2 I 185MR WT N 187MR Compound Hz *2/*4 P 188MR Hz *4 I 189MR Hz *10 I 190MR WT N 191MR Hz *2 I 192MR Hz *4 I 193MR Hz *2 I 194MR Hz *3 I 195MR Hz *17 or *34 I 196MR Hz *4 I 197MR Compound Hz *2/*4 P 198MR Hz *2 I 199MR Compound Hz *2/*4 P ¹P = poor; N = normal; I = intermediate; HI = Higher inducibility; WT = wild type

All of these adolescents were found to have atypical 2D6 haplotypes (see Tables 3-5). The initial results of our analysis of the cytochrome P450 gene alleles in an adult sample also demonstrated a high degree of variability in the polymorphisms of the 2D6 gene (see Tables 3-5). Adolescents with polymorphisms associated with poor 2D6 metabolism were noted to have clinical histories marked by either a poor response to the medication or by a high frequency of side effects when they had been treated with the selective serotonin reuptake inhibitors that are metabolized by the 2D6 enzyme.

Example 2 Determining Genotype of the Dopamine Transporter, Dopamine Receptors, Tryptophan Hydroxylase, Serotonin Receptors, and COMT Genes

The following probes and primers of Table 7 can be generated to detect the genotype of the dopamine transporter (DAT1, SLC6A3), dopamine receptors (DRD1, DRD2, DRD3, DRD4, and DRD5), tryptophan hydroxylase (TPH), serotonin transporter (5-HTT), serotonin receptors (HTR1A, HTR1B, HTR1D, HTR2A, and HTR2C), and COMT genes.

Forward SEQ Symbol Polymorphism Primer Reverse Primer Wt. Probe Variant Probe Stabilizer ID DAT1, >G710A, Q237R CCCCCGCTGAC GACCCCCGAGCCTCACCT AGCTGCCACC AGCTGCCACT GCGGAGGCCCCAG 29-33 SLC6A3 TCCCCTCTG TCC GT >C124T, L42F ACCAAGAGGG CCAATGACGGACAGGAG GAGCTGGTGAG GAGCTGGTGAA CTGCACTCCGTTC 34-38 AAGAAGCACAG AAA TGCTCCT 40 bp VNTR GATGGGGGTC CTGGAGGTCACGGCTCAAG 39-40 CTGGTATGTCT DRD1 >T595G, S199A GTGTGGCATGG GCAATGCGCCGTATTTGT AAAGCTTATTA AGCTTATTACAGC GGATGAGATGGCA 41-45 ACCTTGTCTG TTC CAGA TATGTCCTGC >G150T, R50S CCTGTTTCCTG CCAATACCTGTCCACGCT GGTGTCGGAAC GGTGTCGGAAA CTGATAACGGCAG 46-50 TCGCTGCTC GAT CACAGACC >C110G, T37R GGGCTGGTGGT TCTGACACAGCCAAGGAG CCAGGAGCG CCCAGGAGCC TGGACAGGATGAG 51-55 GGAGAG ATG CAGCGA >A109C, T37P GGGCTGGTGGT TCTGACACAGCCAAGGAG CCAGGAGCGT CCAGGAGCGG GGACAGGATGAGC 56-60 GGAGAG ATG AGCGACA DRD2 >A1051G, TCTCCACAGCA CCGAGAACAATGGCGAG GGTCCGGGT GGTCCGGGC TTTGCCATTGGGC 61-65 T351A CTCCCGACAG CATC ATGGTCTG >C932G, S311C CACCAGCCCA ACCGAGAACAATGGCGA CCATGGTGGG CCATGGTGGC ACGGGTCGGGGA 66-70 CCCGAGAG GCAT GAGTC >C928T, P310S CCGGTACAGC GCCGACTCACCGAGAAC GTGGGACGG GGTGGGACGA GTCGGGGAGAGTC 71-75 CCCATCCC AATG AGCTGG >G460A, V154I GGGGACAGGA GAAGAGGAGTGGGCAGG GGAGATCATGAC TGGAGATCATGAT GGTGACCCGGCGC 76-80 GGGAGGGTTG AGATGGTG TTG AG DRD3 >A25G, S9G AGCAACCAAG GAGCGCGCAGTAGGAGA TTCAGGTGGCT TCAGGTGGCC ACTCAGCTGGCTC 81-85 CCCCAAAGAG GG AGAGATGC DRD4 >T581G, V194G CGCTACAACC GCAGCATGAGCGGGCAG GGACGAGTAGA GGACGAGTAGC CCACGTAGTCGCG 86-90 GGCAGGGTGGG GGTA GTCCTC >C841G, P281A CGCCTTCCCCC GTTGGAGCCGCAGGGGT GGCGCGGG GGCGCGGC GGGCGCACAGTCG 91-95 ACGCCACC CCG GG DRD5 >A889C, T297P TGCGGTTTTCC GGAGTTGGCCCAGCAGAA CCGACAGGGT CGACAGGGG CTTGAGAACCTTG  96-100 ATCAAGAAG GTCTCCTTCTTGAT >G1252A, CAACGCCGAC ACATGCGATCGAAAGGA CGGGGGGAAC CGGGGGGAAT GGCGTTGGGCATC 101-105 V418I TTCTGGAA CTCT ATGTGG >G181A, V61M GCAGTCCAGC GCACGATGGCTGCGGA CTGCGGACAC CTGCGGACAT CAGCACGTTGCCC 106-110 CCGAAATG AG >G185C, C62S CAACGTGCTG CATCTTGCGCTTGTA GGCTGCGC TGGCTGCGG ACACCAGCACGTT 111-115 GTG GCCC >T263G, R88L CGCCAAGATG CGGAATGAAGGAGATGA GCCACGAAGA GCCACGAAGC GGTCTGACACAGC 116-120 ACCAAC CAGAGACA >G1354A, GCAGCTGCCT CATTCGGGGTGAAAGGT GTCCAGCTCC GTCCAGCTCT CAGACAGACTCTG 121-125 W455- ACATCCACAT GTTA CAACAGGGTC HTR1A >G815A, G272D GGAGAGTCGG GAGGCGGGGGCACAAGG GCACAGAGCAC GCACAGAGCAT CCCCAGCCTTGCT 126-130 GGAGCAGGAAC CTCCA >G656T, R219L CTACATCCCGC CCTGCTCCCCGACTCTCCA GCGCAGCTC GCGCAGCTA GGAATATGCGCCC 131-135 TGCTGCTCAT ATAGAGAACCA >C548T, P551L GCTCACTTGGC CATGAGCAGCAGCGGGA CATGCGTCGG CATGCGTCGA GGTCCGAGCGGTC 136-140 TTATTGGCTTC TGT TTCCG CTC >A82G, I28V AGGGCAACAA GCAGCACCAACACCGAC ACGTCGGAGAT CGTCGGAGAC ACCAGTAGTGTTG 141-145 CACCACATCA ACCAT CCGCCG CCAC >G64A, G22S TCAGGGCAAC GCACCAACACCGACACCA GTTGCCGCC GTTGCCGCT GGTCTCAAAGGGA 146-150 AACACCACAT GCCGGT >C47T, P16L TTTTCTTCCCT CCACCACGCACGCATTG AAGGGAGCCG AAAGGGAGCCA GTGGTGATGTGGT 151-155 CCCCCTTCC GTTGTTGCC HTR2A T102C GCTCAACTAC TCACTACGGCTGTCAGTA TTAGCTTCTCCA TTAGCTTCTCCG GAGTTAAAGTCAT 156-160 GAACTCCCTAA AAG TACTGTAGAGCCT GGT HTR1D >C794T, S265L GCTCTGCTCGC GATTTTAGTGGCTTTCCT CCGAGTGCG GCCGAGTGCA AGTGCCCCTCATG 161-165 TCAAC TTC GAGGC HTR2C >C10G, L4V CTCAATTTTAA TTTTGAATAGGAAACACC GCATTCCTCAG GCATTCCTCAC GTTCACCATGATT 166-170 ACTTTGGTTGC CATAA GCTTCAGTCTTAA TTA GC >G68C, C23S ATGATGACAA TCCACCATCGGAGGTATT TCACAGAAATA CACAGAAATATC ATTGCCAAACCAA 171-175 TGATGCTGATG GAAAATG TCAC AG TAGGCCAATTAGG ATGATAATG CYP1A2*IF -164C > A GGCTTCAGTTT TAGGCAAGGGCAGGCAC AAGGAGCTGG GAAGGAGCTGA GTACATGGGGGCC 176-180 CCCGATCTG TAC CCC

Example 3

The genotype of the CYP2D6, CYP2C19, and 5HTTR genes was determined in 97 patients using the methods described above. As indicated in Example 1, CYP2D6*2, *3, *4, *10, *17, and *5 del alleles, CYP2C19*2 (A,B), *3, *4, *5 (A,B), *6, *7, and *8 alleles, and the short/long form of the 5HTTR gene were assessed.

Each patient that was genotyped exhibited either a poor response to at least one antidepressant medication or an intolerable side effect to at least one antidepressant medication. Based on the genotype information, 18 medication profiles were generated (Table 8). In Table 8, PM is poor metabolizer; IM is intermediate metabolizer; EM is extensive metabolizer; s is the short form of the gene; and 1 is the long form of the gene. Each medication profile provides antidepressant medications that are acceptable for use, antidepressant medications to be used with caution, and antidepressant medications to be avoided or used with caution and close monitoring.

TABLE 8 Medication Algorithm Microarray Results from Pharmacological Study Key Gene Genotype Genotype frequencies in Category 2D6 2C19 5HTTTR pharmacological study Med Profile 1 PM PM s/s 2 Med Profile 2 PM PM s/l 2 Med Profile 3 PM PM l/l 1 Med Profile 4 PM EM s/s 1 Med Profile 5 PM EM s/l 6 Med Profile 6 PM EM l/l 1 Med Profile 7 IM PM s/s 1 Med Profile 8 IM PM s/l 7 Med Profile 9 IM PM l/l 3 Med Profile 10 EM PM s/s 0 Med Profile 11 EM PM s/l 6 Med Profile 12 EM PM l/l 3 Med Profile 13 IM EM s/s 2 Med Profile 14 IM EM s/l 17 Med Profile 15 IM EM l/l 16 Med Profile 16 EM EM s/s 7 Med Profile 17 EM EM s/l 10 Med Profile 18 EM EM l/l 12

Medication profile 1 is associated with 2D6 and 2C19 poor metabolizers and the s/s form of the 5HTTR gene. In patients with medication profile 1, bupropion is acceptable for use, mirtazapine and fluvoxamine can be used with caution, while citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, nortriptyline, and sertraline should be avoided or used with caution and close monitoring.

Medication profile 2 is associated with 2D6 and 2C19 poor metabolizers and the s/l form of the 5HTTR gene. In patients with medication profile 2, fluvoxamine and bupropion are acceptable for use, sertraline and mirtazapine can be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

Medication profile 3 is associated with 2D6 and 2C19 poor metabolizers and the 1/1 form of the 5HTTR gene. In patients with medication profile 3, fluvoxamine and bupropion are acceptable for use, sertraline, and mirtazapine can be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

Medication profile 4 is associated with 2D6 poor metabolizers, 2C19 extensive metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 4, bupropion is acceptable for use, mirtazapine, citalopram, fluvoxamine, and escitalopram can be used with caution, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, nortriptyline, and sertraline should be avoided or used with caution and close monitoring.

Medication profile 5 is associated with 2D6 poor metabolizers, 2C19 extensive metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 5, citalopram, fluvoxamine, bupropion and escitalopram are acceptable for use, sertraline and mirtazapine can be used with caution, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

Medication profile 6 is associated with 2D6 poor metabolizers, 2C19 extensive metabolizers, and the l/l form of the 5HTTR gene. In patients with medication profile 6, citalopram, fluvoxamine, bupropion and escitalopram are acceptable for use, sertraline and mirtazapine can be used with caution, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

Medication profile 7 is associated with 2D6 intermediate metabolizers, 2C19 poor metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 7, bupropion is acceptable for use, mirtazapine, venlafaxine, amitriptyline, nortriptyline, fluvoxamine, and sertraline can be used with caution, and citalopram, escitalopram, fluoxetine, imipramine, and paroxetine should be avoided or used with caution and close monitoring.

Medication profile 8 is associated with 2D6 intermediate metabolizers, 2C19 poor metabolizers, and the s/l form of the 5HTTR gene. In patients with medication profile 8, fluvoxamine, bupropion, and sertraline are acceptable for use, mirtazapine, paroxetine, venlafaxine, amitriptyline, and nortriptyline can be used with caution, and citalopram, escitalopram, fluoxetine, and imipramine should be avoided or used with caution and close monitoring.

Medication profile 9 is associated with 2D6 intermediate metabolizers, 2C19 poor metabolizers, and the l/l form of the 5HTTR gene. In patients with medication profile 9, fluvoxamine, bupropion, and sertraline are acceptable for use, mirtazapine, paroxetine, venlafaxine, amitriptyline, and nortriptyline can be used with caution, and citalopram, escitalopram, fluoxetine, and imipramine should be avoided or used with caution and close monitoring.

Medication profile 10 is associated with 2D6 extensive metabolizers, 2019 poor metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 10, bupropion and venlafaxine are acceptable for use, mirtazapine, nortriptyline, imipramine, amitriptyline, sertraline, paroxetine, and fluoxetine can be used with caution, and fluvoxamine, citalopram, and escitalopram should be avoided or used with caution and close monitoring.

Medication profile 11 is associated with 2D6 extensive metabolizers, 2C19 poor met abolizers, and the s/l form of the 5HTTR gene. In patients with medication profile 11, sertraline, venlafaxine, paroxetine, bupropion, and fluoxetine are acceptable for use, and mirtazapine, fluvoxamine, nortriptyline, citalopram, escitalopram, imipramine, and amitriptyline can be used with caution.

Medication profile 12 is associated with 2D6 extensive metabolizers, 2C19 poor metabolizers, and the l/l form of the 5HTTR gene. In patients with medication profile 12, sertraline, venlafaxine, paroxetine, bupropion, and fluoxetine are acceptable for use, and mirtazapine, fluvoxamine, nortriptyline, citalopram, escitalopram, imipramine, and amitriptyline can be used with caution.

Medication profile 13 is associated with 2D6 intermediate metabolizers, 2C19 extensive metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 13, bupropion and mirtazapine are acceptable for use, venlafaxine, amitriptyline, imipramine, nortriptyline, citalopram, fluvoxamine, escitalopram, and sertraline can be used with caution, and paroxetine and fluoxetine should be avoided or used with caution and close monitoring.

Medication profile 14 is associated with 2D6 intermediate metabolizers, 2019 extensive metabolizers, and the s/l form of the 5HTTR gene. In patients with medication profile 14, citalopram, fluvoxamine, bupropion, escitalopram, sertraline, and mirtazapine are acceptable for use, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline can be used with caution.

Medication profile 15 is associated with 2D6 intermediate metabolizers, 2C19 extensive metabolizers, and the l/l form of the 5HTTR gene, hi patients with medication profile 15, citalopram, fluvoxamine, bupropion, escitalopram, sertraline, and mirtazapine are acceptable for use, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline can be used with caution.

Medication profile 16 is associated with 2D6 extensive metabolizers, 2019 extensive metabolizers, and the s/s form of the 5HTTR gene. In patients with medication profile 16, bupropion, mirtazapine, venlafaxine, amitriptyline, imipramine, and nortriptyline are acceptable for use, and paroxetine, fluoxetine, fluvoxamine, sertraline, escitalopram, and citalopram can be used with caution.

Medication profile 17 is associated with 2D6 extensive metabolizers, 2C19 extensive metabolizers, and the s/l form of the 5HTTR gene. Medication profile 18 is associated with 2D6 extensive metabolizers, 2019 extensive metabolizers, and the 1/1 form of the 5HTTR gene. In patients with medication profiles 17 or 18, citalopram, fluvoxamine, bupropion, escitalopram, sertraline, mirtazapine, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline are acceptable for use.

Example 4

The genotype of the CYP2D6, CYP2C19, CYP3A4, CYP1A2, and HTR2A genes was determined in ten patients using the methods described above. The genotype of the 5HTTR gene was determined by amplifying the polymorphic region of the gene using primers flanking the region. Reaction products were sized using electrophoresis (Agilent Technologies). CYP3A4*1B, *2, *5, *6, *12, *13, *15A, *17, and *18A alleles, CYP1A2*1 F allele, and the MICA 102 polymorphism were detected

Based on the CYP2D6 genotype, patients were classified as being poor, intermediate, extensive, or ultra rapid metabolizers. For each of the CYP2C19, CYP3A4, and CYP1A2 genotypes, patients were classified as being either extensive or not extensive (i.e., poor) metabolizers. For the 5HTTR genotype, patients were classified as having the short/short (s/s) form of the gene, the short/long (s/l) form of the gene, or the long/long (l/l) form of the gene. For the HR2A genotype, patients were classified as having either the C/C allele or as having the T/T or T/C alleles.

The following six rules were used to sort the genotypes into 192 treatment recommendations. The first rule relates to the CYP2D6 genotype. For poor metabolizers and ultra rapid metabolizers, CYP2D6 substrates were placed in the “use very cautiously with close monitoring” category. For intermediate metabolizers, exclusive CYP2D6 substrates were placed in the “avoid or use very cautiously with close monitoring” category. For intermediate metabolizers, substrates that are partially metabolized by CYP2D6 were placed in the “use cautiously” category. For extensive metabolizers, CYP2D6 substrates were placed in the “acceptable for use” category.

After applying the first rule, the second rule relating to the 2C19 genotype was applied. For poor metabolizers based on their CYP2C19 genotype, exclusive CYP2C19 substrates were placed in the “avoid or use very cautiously with close monitoring” category. For extensive or intermediate metabolizers based on their CYP2C19 genotype, CYP2C19 substrates were placed in the “acceptable for use” category. Substrates that are primarily metabolized by both CYP2D6 and CYP2C19 were placed in the “avoid or use very cautiously with close monitoring” category if CYP2C19 was poor and CYP2D6 was either intermediate or poor, or placed in the “use cautiously” category if CYP2C19 was poor and CYP2D6 was extensive.

A third rule relating to the 3A4 genotype was applied after the first and second rules were applied. For poor metabolizers based on their CYP3A4 genotype, exclusive CYP3A4 substrates were placed in the “avoid or use very cautiously with close monitoring” category. For extensive or intermediate metabolizers based on their CYP3A4 genotype, CYP3A4 substrates were placed in the “acceptable for use” category. Substrates that are primarily metabolized by any two of the following three enzymes: CYP2D6, CYP2C19, or CYP3A4, were placed in the “avoid or use very cautiously with close monitoring” category if any two of the three were intermediate or poor metabolizers or placed in the “use cautiously” category if one of the two most primary enzymes was a poor metabolizer.

The fourth rule relating to the 1A2 genotype was applied after applying the first three rules. For poor metabolizers based on CYP1A2 genotype, CYP1A2 substrates were placed in the “avoid or use very cautiously with close monitoring” category (include fluvoxamine as a primary CYP1A2 substrate). For poor metabolizers based on the CYP1A2 genotype, substrates having the CYP1A2 pathway as one of their primary pathways were placed in the “use with caution” category.

After applying the first four rules, the fifth rule relating to the 5HTTR genotype was applied. For the homozygous short genotype (s/s) of the serotonin transporter gene, SSRI's that were categorized as “acceptable for use” were shifted to the “use with caution” category. For the homozygous short genotype (s/s) of the serotonin transporter gene, SSRI's that were categorized as “use with caution” were shifted to the “avoid or use very cautiously with close monitoring” category. Antidepressants, outside the SSRI class, do not shift category based on 5HTTR genotype.

Rule six, which relates to the HR2A genotype, was applied after applying the first five rules. For the “mutant form” of the 5-hydroxytryptamine serotonin receptor 2A gene, if paroxetine was in the “acceptable for use” category, it was shifted to the “use with caution” category. If paroxetine was in the “use with caution” category, it was shifted to the “avoid or use very cautiously with close monitoring.

Based on the genotype of the patients and associated medication profile, medications appropriate for the patient were selected. Table 9 summarizes the data for the 10 patients. Case studies of four of the seven patients are presented below.

The first case is a 16-year-old white female who was genotyped and found to be a poor metabolizer for CYP2D6 (*3/*9) as well as being a poor metabolizer for CYP2C19 (*2/*2). She was diagnosed at 6 years of age as having ADHD and treated unsuccessfully. She subsequently developed symptoms of depression and mood dysregulation and was placed on many medications that were metabolized by CYP2D6, such as EFFEXOR. Additionally, she was homozygous for a short form of the serotonin transporter gene. Based on her specific genotype of these key genes, the algorithm demonstrates that there are a relatively small number of potential strategies for controlling her depression. The medication profile generated based on her genotype of HR2A TT/TC, s/s form of 5HTTR, extensive metabolizer for 1A2 and 3A4, and poor metabolizer for 2C19 and 2D6 indicated that bupropion was acceptable for use, mirtazapine and fluoxamine should be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, nortriptyline, and sertraline should be avoided or used with caution and monitoring.

Before genotyping, she had a particularly unusual side effect to citalopram in that she developed hallucinations (incidence less than 1%). This atypical and dangerous side effect could have been avoided if she had been genotyped prior to receiving this medication. The metabolic enzyme that is primarily involved in the metabolism of citalopram is CYP2C19 and a secondary metabolic pathway is CYP2D6. Both of these pathways were nonfunctional in this patient. Therefore, even on a modest dose of citalopram, this patient developed acute auditory and visual hallucinations, which were subsequently treated with an atypical antipsycliotic medication. Additionally, based on an inability to develop an effective treatment plan for this patient, she required three acute hospitalizations and a partial hospitalization experience prior to being referred to a state hospital for two months. She was unsuccessfully tried on more than a dozen psychotropic medications. There is no question that her care could have been greatly improved and perhaps some of her hospitalizations could have been avoided if clinicians had been aware of her inability to use many of the drugs that had been prescribed at their standard doses. Additionally, her serotonin transporter genotype (C/C) would have suggested a poor response to many of the medications that had been prescribed. Subsequent to genotyping, she was placed on a combination of two mood stabilizers that are metabolized by an alternative pathway (LAMICTAL (lamotrigine) and TRILEPTAL (oxcarbazepine)). She is also tolerating treatment with ABILIFY (aripiprazole), which in addition to being metabolized by CYP2D6 is metabolized by CYP3A4, for which she has an internal metabolic pathway. She is now doing reasonably well.

The second case is a 42-year-old Caucasian female with a persistent depressive illness and pervasive anxiety and a long-standing history of psychotropic medication intolerance. This patient had a comprehensive neurological evaluation, which failed to provide an organic explanation for her symptoms. This patient has been treated with at least nine psychotropic medications in the past. Striking examples of CYP2D6 intolerance include the development of acute hallucinations while taking paroxetine and a very low tolerance to amitriptyline. This patient had poor CYP2D6 metabolism and intermediate CYP2C19 metabolism. This represents a rare genotype that would suggest difficulties in tolerating treatment with citalopram. Indeed, she did have irritability and lethargy when taking citalopram. To compound this problem, she also was homozygous for the short form of the serotonin transporter, which is associated with a poor response to serotonin receptive inhibitors. She had experienced multiple outpatient visits for her psychiatric problems. Awareness of her genotype earlier in her history would clearly have allowed her clinician to avoid a number of medications that she could not tolerate. The medication profile generated based on her genotype of HR2A TT/TC, s/s form of 5HTTR, extensive metabolizer for 1A2 and 3A4, and poor metabolizer for 2C19 and 2D6 indicated that bupropion was acceptable for use, mirtazapine and fluoxamine should be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, nortriptyline, and sertraline should be avoided or used with caution and monitoring.

The third case history is a 67-year-old woman with a diagnosis of major depressive disorder, who had been treated with seven different medications in the past. These had not been effective, and she had required two psychiatric hospitalizations. It was determined that she was a poor metabolizer for both CYP2D6 (*4/*9) and CYP2C19 (*2/*2). Additionally, she was homozygous for the long form of the serotonin transporter gene. Based on her genotype of HR2A C/C, l/l form of 5HTTR, extensive metabolizer for 1A2 and 3A4, and poor metabolizer for 2C19 and 2D6, the associated medication profile indicated that fluoxamine and bupropion were acceptable for use, sertraline and mirtazapine should be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and monitoring. If her genotype had been known, it would have been possible to anticipate that she would have had side effects to both EFFEXOR (Venlafaxine) and LEXAPRO (escitalopram oxalate). This could have been predicted, given that two of the three metabolic pathways involved in the metabolism of LEXAPRO were completely inactive. Subsequent to her genotyping, she was successfully treated with electroconvulsive therapy (ECT) and concurrently successfully treated with REMERON (mirtazapine), trazodone, and KLONOPIN (clonazepam), which have alternative pathways of metabolism to CYP2D6 and CYP2C19.

The fourth case is a 50-year-old Caucasian woman who had been diagnosed as having dysthymia and chronic anxiety. She had been placed on multiple psychotropic drugs, including two CYP2D6 medications that she could not tolerate. Her genotype indicated that, she was a poor CYP2D6 metabolizer (*4/*4) and an intermediate CYP2C19 metabolizer (*1/*2). She was heterozygous for the long and short form of the serotonin transporter gene. Based on her genotype of ER2A C/C, l/s form of 5HTTR, extensive metabolizer for 1A2 and 3A4, and poor metabolizer for 2C19 and 2D6, the associated medication profile indicated that fluoxamine and bupropion were acceptable for use, sertraline and mirtazapine should be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and monitoring.

This patient had multiple hospitalizations at considerable cost. She also found it difficult to tolerate a range of nonpsychiatric medications, including being able to take only one dose of PROZAC (fluoxetine) and ZOLOFT (sertraline) before finding that she was acutely uncomfortable. Of particular interest, she ultimately responded to a subclinical dose of paroxetine, which represents an example of using the algorithm to find a method to administer an effective medication that in normal doses could not have been tolerated.

TABLE 9 Sub. HTR2A No Age Gender Race Diagnosis 2D6 2C19 3A4 1A2 51HTTR T102C 1 16 F Cauc. ADHD; adjustment 3/9 Poor 2/2 Poor *1*1 *1/*1F s/s TC disorder NOS; Depression Freq: Freq: NOS; Psychosis NOS; *3 = 1% *2 = Bipolar Disorder *9 = 3% 14% 2 42 F Cauc. Generalized anxiety 4/4 Poor 1/2 Poor *1*1 *1/*1F s/s TC disorder, Depressive Freq: Freq: Disorder NOS; Panic *4 = 18% *1 = 86% disorder *2 = 14% 3 67 F Other Major depressive 4/9 Poor 2/2 Poor *1*1 *1/*1 l/l CC disorder Freq: Freq: *9 = 3% *2 = 14% 4 50 F Cauc. Panic disorder; 4/4 Poor 1/2 Int *1*1 *1/*1F s/l CC generalized anxiety Freq. Freq: disorder; dysthymia, *4 = 18% *1 = 86% somatogorm disorder *2 = 14% NOS 5 39 F Cauc. Generalized anxiety 3/4 Poor 1/2 Poor *1*1 *1/*1F l/l CC disorder Freq. Freq: *3 = 33% *1 = 86% *4 = 18% *2 = 14% 6 12 F Hispanic Bipolar disorder NOS 1/1 ext. 1/2 Int. *1/*1B *1/*1 s/l CC freq. Freq. *1 = 37% *1 = 86% *2 = 14% 7 54 M Bipolar 2P/2P ext. 1/1 Ext. *1*1 *1/*1F s/l TC disorder freq: Freq: *2P = unk. *1 = 86% 8 47 F Cauc. Major depressive 1/5 Int. 1/1 Ext. *1/*1B *1/*I l/l TC psychotic *1 = 37% *1 = 37% *5 = 4% 9 25 F Canc. Panic disorder 2/2P ext. 1/1 Ext. *1*1 *1/*1 s/1 TC freq: Freq: *2 = 33% *1 = 37% *2P = unk. 10 16 F Asian- Major depressive disorder 4/10 Poor 1/1 Ext. *1*1 *1F/*1F s/l TT American vs. bipolar Freq: Freq: *4 = 18% *1 = 86% *10 = 2%

Based on patient no. 5's genotype of HR2A C/C, 1/1 form of 5HTTR, extensive inetabolizer for 1A2 and 3A4, and poor metabolizer for 2C19 and 2D6, the associated medication profile indicated that fluoxamine and bupropion were acceptable for use, sertraline and mirtazapine should be used with caution, and citalopram, escitalopram, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and monitoring.

The medication profile associated with patient no, 6's genotype (HTR2A CC; 5HTTR l/s; 1A2, 3A4, and 2D6 extensive; and 2C19 poor) indicates that sertraline, venlafaxine, bupropion, and fluoxetine are acceptable for use, and mirtazapine, fluvoxamine, nortriptyline, citalopram, escitalopram, imipramine, amitriptyline, and paroxetine can be used with caution.

The medication profile associated with patient no, 7's genotype (HTR2A TT/TC; 5HTTR l/s; 1A2, 3A4, and 2019 extensive; and 2D6 ultra rapid) indicates that citalopram, fluvoxamine, bupropion, and escitalopram are acceptable for use; sertraline and mirtazapine can be used with caution, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

The medication profile associated with patient no. 8's genotype (HTR2A TT/TC; 5HTTR l/l; 1A2, 3A4, and 2C19 extensive; and 2D6 intermediate) indicates that citalopram, fluvoxamine, bupropion, escitalopram, sertraline, and mirtazapine are acceptable for use and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

The medication profile associated with patient no. 9's genotype (HTR2A TT/TC; 5HTTR l/s; and 1A2, 3A4, 2C19, and 2D6 extensive) indicates that citalopram, fluvoxamine, bupropion, escitalopram, sertraline, mirtazapine, paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline are acceptable for use.

The medication profile associated with patient no, 10's genotype (HTR2A TT/TC; 5HTTR l/s; 1A2, 3A4, and 2C19 extensive; and 2D6 poor) indicates that citalopram, fluvoxamine, bupropion, and escitalopram are acceptable for use, sertraline and mirtazapine should be used with caution, and paroxetine, fluoxetine, venlafaxine, amitriptyline, imipramine, and nortriptyline should be avoided or used with caution and close monitoring.

OTHER EMBODIMENTS

It is to be understood that while the invention has been described in conjunction with the detailed description thereof, the foregoing description is intended to illustrate and not limit the scope of the invention, which is defined by the scope of the appended claims. Other aspects, advantages, and modifications are within the scope of the following claims 

1-47. (canceled)
 48. A computer-implemented method for selecting an antidepressant medication for a patient in need thereof, the method comprising identifying the patient's clinical phenotype as depression based on data received by the computer system, receiving a set of medications used to treat depression, associating each of the medications with a combination genotype and a category of use to create at least one medication profile for each medication, the combination genotype comprising the genotypes of a panel of at least three genes, a cytochrome P450 CYP2D6 gene, a cytochrome P450 CYP2C19 gene, and a serotonin transporter gene 5HTTR, storing the medication profiles and genotypes such that a medication profile associated with a specific genotype can be identified, receiving the patient's combination genotype for the panel of at least three genes, identifying from the stored medication profiles at least one identified medication profile associated with the patient's combination genotype, and outputting a patient-specific report comprising the at least one identified medication profile associated with the patient's genotype.
 49. The method of claim 48, wherein the medications are selected from citalopram, fluvoxamine, bupropion, escitalopram, sertraline, mirtazapine, fluoxetine, venlafaxine, amitriptyline, imipramine, paroxetine, and nortriptyline.
 50. The method of claim 49, wherein the panel of at least three genes comprises (i) a cytochrome P450 CYP2D6 gene consisting of two alleles independently selected from the group consisting of fully functional alleles 2D6*1A, *2, *2N, and *35, partially functional alleles 2D6*10 and *17, and non-functional alleles 2D6*3, *4, *5, *6, *7, *8, and *12, (ii) a cytochrome P450 CYP2C19 gene consisting of two alleles independently selected from the group consisting of fully functional alleles 2C19*1A and *1B, and non-functional alleles 2C19*2A, *2B, *3, *4, *5A, *5B, *6, *7 and *8, and (iii) a serotonin transporter gene 5HTTR consisting of two alleles independently selected from the short (s) and long (1) forms of the 5HTTR gene.
 51. The method of claim 50, further comprising assigning a metabolic phenotype to each cytochrome P450 gene according to the number of functional alleles in the genotype, the metabolic phenotype selected from poor (P), intermediate (I), and extensive (E).
 52. The method of claim 51, wherein the at least one medication profile has a data structure comprising<medication::2D6/2C19 metabolic phenotype; 5HTTR genotype::use category>and is selected from the group consisting of: amitriptyline::P/P or P/E; s/s, s/l or l/l::avoid or use with caution and close monitoring amitriptyline::I/P, E/P or I/E; s/s, s/l or l/l::use with caution amitriptyline::E/E; s/s, s/l or l/l::acceptable for use as directed; bupropion::P/P, P/E, I/P, E/P, I/E, or E/E; s/s, s/l or l/l::acceptable for use as directed; citalopram::P/P or I/P; s/s, s/l or l/l::avoid or use with caution and close monitoring citalopram::E/P; s/s::avoid or use with caution and close monitoring citalopram::P/E, I/E or E/E; s/s::use with caution citalopram::E/P; s/l or l/l::use with caution citalopram::P/E, I/E or E/E; s/l or l/l::acceptable for use as directed; escitalopram::P/P or I/P; s/s, s/l or l/l::avoid or use with caution and close monitoring escitalopram::E/P; s/s::avoid or use with caution and close monitoring escitalopram::P/E, I/E or E/E; s/s::use with caution escitalopram::E/P; s/l or l/l::use with caution escitalopram::P/E, I/E or E/E; s/l or l/l::acceptable for use as directed; fluoxetine::P/P, P/E or I/P; s/s, s/l or l/l::avoid or use with caution and close monitoring fluoxetine::I/E; s/s::avoid or use with caution and close monitoring fluoxetine::E/P or E/E; s/s::use with caution fluoxetine::I/E; s/l or l/l::use with caution fluoxetine::E/P or E/E; s/l or l/l::acceptable for use as directed; fluvoxamine::E/P; s/s::avoid or use with caution and close monitoring fluvoxamine::P/P, P/E, I/P, I/E or E/E; s/s::use with caution fluvoxamine::E/P; s/l or l/l::use with caution fluvoxamine::P/P, P/E, I/P, I/E or E/E; s/l or l/l::acceptable for use as directed; imipramine::P/P, P/E or I/P; s/s, s/l or l/l::avoid or use with caution and close monitoring imipramine::E/P or I/E; s/s, s/l or l/l::use with caution imipramine::E/E; s/s, s/l or l/l::acceptable for use as directed; mirtazapine::P/P, P/E, I/P or E/P; s/s, s/l or l/l::use with caution mirtazapine::I/E or E/E; s/s, s/l or l/l::acceptable for use as directed; nortriptyline::P/P or P/E; s/s, s/l or l/l::avoid or use with caution and close monitoring nortriptyline::I/P, E/P or I/E; s/s, s/l or l/l::use with caution nortriptyline::E/E; s/s, s/l or l/l::acceptable for use as directed; paroxetine::P/P or P/E; s/s, s/l or l/l::avoid or use with caution and close monitoring paroxetine::I/P or I/E; s/s::avoid or use with caution and close monitoring paroxetine::I/P or I/E; s/l or l/l::use with caution paroxetine::E/P or E/E; s/s::use with caution paroxetine::E/P or E/E; s/l or l/l::acceptable for use as directed; sertraline::P/P or P/E; s/s::avoid or use with caution and close monitoring sertraline::P/P or P/E; s/l or l/l::use with caution sertraline::I/P, E/P, I/E or E/E; s/s::use with caution sertraline::I/P, E/P, I/E or E/E; s/l or l/l::acceptable for use as directed; venlafaxine::P/P or P/E; s/s, s/l or l/l::avoid or use with caution and close monitoring venlafaxine::I/P or I/E; s/s, s/l or l/l::use with caution venlafaxine::E/P or E/E; s/s, s/l or l/l::acceptable for use as directed;


53. The method of claim 48, wherein the category of use is selected from acceptable for use as directed, use with caution, avoid or use with caution and close monitoring.
 54. The method of claim 48, further comprising querying a database of genotypes.
 55. The method of claim 48, further comprising querying a database of electronic health records.
 56. The method of claim 48, further comprising synthesizing the genotypes, phenotypes, and medication profiles into an individualized composite report.
 57. The method of claim 48, wherein the patient's genotype is received directly from equipment used in determining the patient's genotype.
 58. The method of claim 48, wherein the medication profile for the patient includes a ranking of several medications.
 59. The method of claim 48, wherein the ranking is based on patient-specific co-factors. 